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1. Measurement Scheme 

In using mass spectrometers to measure abun- 
dance ratios of isotopes in an element, a general 
procedure is to relate all minor abundant isotopes 
to one major abundant isotope at the same point in 
time. Ion currents in a scan of the isotopic ion 
beams are either adjusted for time differences or 
normalized by a symmetrical arrangement of mea- 
surement sequence followed by averaging. 

This paper proposes a scheme of measurement 
sequences that has been used for the redetermina- 
tion of the molar mass (atomic weight) of silicon at 
the Central Bureau for Nuclear Measurements 
(CBNM). 1 This scheme avoids correlations among 



1 The molar mass of a particular specimen of silicon thus deter- 
mined was used in a project, in collaboration with the 
Physikalisch-Technische Bundesanstalt (PTB), Braunschweig, 
Germany, for the redetermination of the Avogadro Constant 
[1,2]. 



the measured ratios caused by normalizing all ion 
current measurements to that of the largest ion 
current. It also provides additional information for 
checking on the consistency of these ratios within 
a cycle of scans. Measurements of isotope abun- 
dance ratios of silicon are used below as an illustra- 
tion. 

Recalling that there are three stable isotopes of 
mass numbers 28, 29, and 30, with atomic masses 
M ( 28 Si), M( 29 Si), and M ( 30 Si) known almost exactly 
[3], the molar mass of a sample of silicon, M, can 
be determined as: 

M(Si) = XM(Si)f i 

where fi is the fractional isotope abundance of iso- 
tope 'Si, with 2/i = l. 
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In our experiment, measurements were made on 
SiF4 gas samples by means of a gas mass spectrom- 
eter equipped with a Faraday-cup collector. Ion 
currents were measured for 28 SiF3 + (mass position 
85), 29 SiF 3 + (86), and 30 SiF 3 + (87). Ion current ra- 
tios, R(i/j), are formed from these measurements. 
These ratios, after corrections by proportionality 
factors, are estimates of ratios of fractional abun- 
dances,//^. Denoting the corrected ratios as r(i/j) 
and equating these to the fi/fj, it follows that 



/28-j 



+ r (29/28) +r (30/28)' 



(1) 



and 



f»=r(29J28)fa, 
f 3 o=r (30/28)/*. 

Scaling the conversion of ion-current ratios into 
isotope-abundance ratios, i.e., determining the pro- 
portionality factors is accomplished by means of ac- 
curately synthesized isotope mixtures with isotope 
abundance ratios close to those of the unknown 
silicon samples. 

In previous work [1], a two-ratio scheme was 
used, because Eq. (1) requires only the measured 
values of r(29/28) and r(30/28), i.e., the ratios of 
minor abundant isotopes (3 + to 4 + %) to the major 
abundant isotope 28 Si (92 + %) The sequence of 
measurements in a scan proceeds as depicted in 
Fig. 1. There is one average for each intensity mea- 
sured: 7(1) for 28 SiF 3 + , 7(2) for 29 SiF 3 + , and 7(3) for 
3U SiF 3 + . Thus the ratios are: 



7?(86/85) =7(2)/I(l); 7?(87/85) =7(3)/7(l). 

Eight to ten scans are measured during one cycle. 
A linear least squares fit of these measured ratios 
vs time would yield estimates of the ratios at t , the 
instant sample gas is introduced into the ionizing 
chamber, to calculate and to allow for mass frac- 
tionation effect. The residual standard deviation, 
and the relative standard deviations of the 7?'s are 
also computed for each cycle. 

The two-ratio scheme has been in use and ap- 
pears to work well. However, two questions remain 
unanswered: 

1. Since both 7? (86/85) and 7? (87/85) used the 
same value of 7(85) in the denominator, the 
two ratios are correlated. What is the effect of 
correlation on the resulting calculated / 2 s, /», 
and /30? In addition, if the intensity of the ma- 
jor isotope is measured higher (or lower) than 
it should, because of minor nonlinearity in the 
instrumentation system, then the correspond- 
ing abundance will also be higher (or lower). 

2. How do we know if the intensities measured in 
a scan, or in a cycle, are correct or consistent? 

To answer these questions, a three-ratio scheme is 
devised such that all three ratios, 7? (87/85), 7? (85/ 
86), and 7? (86/87), are measured, using six inde- 
pendently measured averages of intensities. Thus 
12 ion intensities are measured in a sequence in 
one scan as shown in Fig. 2. For the scheme shown, 
there are two average intensities for 85, 7(1) and 
7(4), two averages for 86, 7(3), and 7(6), 



Mass Position 



1(3) 



1(2) 

— t— 



R(86/85) = I(2)/I(1) 
R(S7/85) = I(3)/I(1) 



10) 



Fig. 1. Two ratio scheme. One scan of a two ratio scheme is shown 
here where two measurements of ion currents of eaeh isotopic ion 
beam are made, centered about /. Averages of those three pairs of 
measurements, 1(1), 1(2), and 7(3), are used to form the two ratios 
R (86/85) and R (87/85). 
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Mass Position 85 87 



85 87 



86 87 85 86 87 85 



1(6) 

■H— 

K5) 



R(87/85) = I(2)/1(1) 
R(85/86) = I(4)/I(3) 
R(86/87) = I(6)/I(5) 



1(4) 

-)— 

K3) 



K2) 

— \— 



10) 

Fig. 2. Three ratio scheme. One scan of a three ratio scheme is shown 
in this diagram where four measurements of ion currents of each ion 
beam are made, all centered about i. Averages of each two measure- 
ments equidistant about /, 1(1) through 1(6), are used to form the three 
ratios R (87/85), R (85/86) and R (86/87). 



and two averages for 87, 1(2), and 1(5). The ratios 
may be formed as 

i?(87/85)=/(2)/I(l), R (85/86) =/(4)/I(3), and 

R (86/87) =/(6)/I(5). 

Hence, there is no obvious correlation among the 
ratios to worry about because all ratios are formed 
with independently measured intensities. Further- 
more, the "redundant" third ratio R (86/87) can be 
used to form a "closure" to check on the consis- 
tency of the measured ratios [4], i.e., 

R (87/85) R (85/86) R (86/87) = 1, 

if there are no measurement errors. With measure- 
ment errors, we may formulate 

R (87/85) R (85/86) R (86/87) = 1 + e, 

and use the e values as a control on the measured 
ratios. We note that 



Average e=0 



Standard 

i?'s) 2 } !/2 . 



deviation of e~ {^(relative s.d. of 



The relative standard deviations, computed from 
the residual standard deviations resulting from the 
linear fit of R's in a cycle, can be used in the above 
expression. 

A control chart on e can then be constructed to 
monitor measurements within a run. If e's are pre- 
dominately positive (negative), it is an indication of 
the presence of a systematic error at some, identifi- 
able, point in the measurement procedure. Investi- 
gation as to its cause is in order. 

The three-ratio scheme has been implemented in 
this laboratory (CBNM) for about 6 months now 
and seems to work well [2]. With computer con- 
trolled measurement and summary, the additional 
work is minimal once the software is prepared. A 
typical data summary sheet for the measurement of 
a silicon specimen is shown in Table 1, listing the 
three ratios measured in a cycle of ten scans, to- 
gether with the e's calculated for each scan, extrap- 
olated values at t and the relative standard 
deviations. 
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Table 1. Typical data summary sheet for a cycle of ten scans. Observed ratios of ion currents R(86/85), #(85/87) and 
R (87/86) arc listed in columns 2, 3, and 4, respectively, for the ten scans. The (1 + e) values area calculated for each scan and 
for the extrapolated values at to 



Scan no. 



/( 2S Si)//( 2B Si) 



/(^Siy/fsi) 



/( 3 °Si)//( 2 "Si) 



1+<T 



1 

2 
3 
4 
5 
6 
7 
8 
9 
10 



0.050 5975 


29.8800 


0.661 443 


1.000 005 


0.0505971 


29.871 2 


0.661 547 


0.999 860 


0.050 606 


29.864 9 


0.661 637 


0.999 961 


0.050 615 4 


29.858 5 


0.661 706 


1.000 036 


0.050 6175 


29.849 7 


0.661 798 


0.999 922 


0.050 622 9 


29.842 2 


0.661 888 


0.999 913 


0.050 628 5 


29.834 4 


0.661 954 


0.999 862 


0.050 636 2 


29.826 7 


0.662 056 


0.999 910 


0.050 642 7 


29.819 4 


0.662 123 


0.999 895 


0.050 645 7 


29.812 9 


0.662 238 


0.999 910 



Mean 
Std. dev. 






0.999 927 
57 


Extrapolation 


Value (<=/o) 
Std. dev. 


0.0505745 
20 


29.906 7 0.661 149 
7 11 


0.999 998 
48 


Rel. std. dev. 


3.9 -lO" 5 


2.3 -lO -5 1.7 -10" 5 


4.8 . 10" 5 



With the addition of the third ratio, there are 
also a number of features that can be used to check 
on the accuracy of the mass spectrometric mea- 
surements: 

a. We note that in Fig. 2, the three ratios can be 
computed by different pairing of the intensities, 
e.g., 

J? (87/85) =/(2)//(4), R (85/86) = /(l)//(6), 
R (86/87) =7(3)// (5), 

and so on. Ratios computed from different pairings 
can be compared. If they differ consistently, the 
cause should be investigated. 

b. If we denote 85 by A, 86 by B and 87 by C, 
then the sequence in Fig. 2 can be represented as 

ACBACBBCABCA. 

By permuting the positions of these isotopes, we 
could use also 

ABCABCCBACBA 
BACBACCABCAB 
BCABCAACBACB 
CBACBAABCABC 
CABCABBACBAC. 

The essential difference of these six sequences is 
the position of the major isotope A relating to the 
minor isotopes. If these six sequences yield ratios 
that are different beyond experimental errors, it is 



an indication that adjustments should be made on 
mass position, interference, or other factors. 

c. The ratio of the minor isotopes, R (86/87) is 
much more sensitive to small changes than the 
other two ratios where the major peak dominates 
the behavior of these ratios. For example, when a 
natural silicon sample is measured after an en- 
riched 29 Si, the mass spectrometer seems to remem- 
ber the last measurement (by adsorption to the 
walls or other reasons) and yields a higher 86 inten- 
sity than actually present in the natural silicon. This 
effect shows more clearly in the ratio (86/87) than 
the other two. Hence it can be used to check 
whether there has been enough flushing and clean- 
ing of the ion source to erase the memory. 

d. If we denote the three measured ratios byR u 
i?2, and i?3, then the least squares adjusted ratios are 
[5]: 

ri -\RiRj ;r2 to; ;r3 tei • 

ri and r 2 can be used in Eq. (1). The relative stan- 
dard deviations of the r's is about 0.82 (V2/3) of 
those for the R 's. 

e. If a mass spectrometer has three Faraday 
cups to measure the three intensities separately (as 
is the case of the new spectrometer at CBNM), a 
similar scheme may be devised to yield three inde- 
pendent ratios as shown in Fig. 3. For each isotope, 
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Cup 1 



Cup 2 



Cup 3 



_J I 1 1 1 — p. 

a b | b a 

-J 1 1 1 1 — ^ 

a bib a 



■>85 
b85 



R(87/85) = a 87 / b 85 
R(85/86) = S S5 / b 86 
R(86/87) = a 86 / b 87 



ajv 
b 87 



Fig. 3. Scheme for three Faraday cups. One scan of a scheme using 
three Faraday cups is shown in this diagram. Two averages of ion cur- 
rents of each ion beam, a and b, are used to form the three ratios 
tf (87/85), R (85/86), and tf (86/87). 



e.g., 85 using cup 1, four ion currents are measured 
at equal intervals in time in a scan. The two end 
ones are designated by a 's, and the two middle ones 
by b's. Similarity, for 86 using cup 2, and 87 using 
cup 3, measured during the same time intervals as 
those of 85. Thus for each isotope, we have two 
averages, a and b , both corresponding to time 1. 
Hence, three independent^ratios can be computed, 
using the six averages ass, oss, a~s6, 086, ~a&, and bun, as 
follows: 

R (86/85) =£* ; R (85/87) =|^ ; R (87/86) =§E , 

085 087 086 

Here the e's calculated from the three independent 
ratios would also check the consistency of the three 
cups. 

With the objective to determine the molar mass 
of silicon accurate to one part in 10 7 , thus requiring 
a precision in ratio measurements to parts in the 10 5 
range, it is imperative to investigate all avenues of 
improvements. The three-ratio scheme provides 
symmetry and redundancy and appears to be a 
helpful step in this direction. 

In the above, we have used the three isotopes of 
silicon as an example to illustrate the three-ratio 
scheme. It is obvious that any three- or more- iso- 
topes of a polyisotopic element can be treated in 
the same manner. The selection of the particular set 
of three isotopes is, of course, a decision the exper- 
imenter must make to suit his objectives. 
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